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> 1 I. INTRODUCTION 

in 

I Exclusion statistics is one possible way to generalize the pattern of bosonic or fermionic 

q ■ particles. A creation operator a> can create at most k — 1 particles at a given site (or 

mode) by stating fc-nilpotency: acting on a vacuum |0), 

( a t) n |0) ^ Ofarn = 0,--,jfc-l 
.W. (« t ) fc |0> = 0- W 

Standard bosons are recovered in the limit k — > oo while the fermionic Pauli exclusion 
principle corresponds to k — 2. 

Amongst proposals of systems satisfying exclusion statistics with finite k > 3, we are 
interested in this work in the so called fc-nilpotent g-boson particles^. We first review 
the operator formulation of quantum mechanics for one degree of freedom and for a 
system of m independent g-boson modes. Then we discuss the construction of coherent 
states, for which it is necessary to introduce /c-nilpotent para-Grassmann numbers (cf. 
Grassmann numbers in the fermionic case). The order of nilpotency k is usually related 
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to commutation rules for para-Grassmann numbers^ (in a similar way that fermionic 
k = 2 nilpotency leads to anticommuting Grassmann numbers). Ordering issues are 
then important and have made cumbersome the manipulation of coherent states. Several 
authors have dealt with this problem^- generating a variety of conventions and finding 
unavoidable difficulties, in particular when dealing with multiparticle states. 

We explore in this work the consistency of using non-standard para-Grassmann commu- 
tation relations and a normal order convention^, highly simplifying the operations but 
still retaining the essence of exclusion statistics. We then write down a para-Grassmann 
symbolic expression for the trace of operators and use it to study the thermodynamics 
of g-boson systems. The partition function for non-interacting systems can be readily 
computed, alongside with derived quantities like the mean free energy and the specific 
heat, observing that fc-nilpotent bosons interpolate the features of Fermi-Dirac and Bose- 
Einstein statistics. 

The present formalism could be useful in the study of strongly correlated systems in low 
dimensions, where effective quasi-particles with exclusion statistical properties seem to 
be ubiquitous. The constraints on available states for spin S particles, or for electrons in 
t-J models, or for fermionic and bosonic occupation in representations of spin operators, 
impose rules on statistical distributions which often manifest in fractional statistics. No- 
ticeably, the "fractional exclusion" statistics characterized by Haldane^ is realized in sev- 
eral strongly correlated systems in one and two dimensions like the Haldane-Shastry spin 
chain^ and generalizations^ and the Fractional Quantum Hall Effect^ (where fractional 
exclusion statistics is consistent with anyon braiding statistics^). In one dimensional 
Conformal Field Theories, the underlying Yangian symmetry allows for the construction 
of a basis of quasi-particle excitations which also have been proved to obey of exclusion 
statistics^. 

II. Q-BOSON OPERATORS 

The origin of g-bosons finds its roots in a Schwinger-like bosonic representation of the 
quantum deformed SU(2) q generators^, where q is a real deformation parameter; later, 
the consideration of q as a rational phased led to nilpotent operators. 

We consider in this work a set of m, fc-nilpotent, independent g-boson modes a^, aj 
(i = 1, ■ ■ • ,m). For each mode^ aj is an annihilation operator and a], its Hermitian 
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conjugate, is the corresponding creation operator satisfying the q- commutation relations 

aid\ — qa\ci{ = q~ Ni , (2) 

and conjugate relations 

did] — g _1 a|aj = q Ni , (3) 

where q = e 1 ^ , k G N, k > 2, and Ni is a number operator that, from (j2J), can be related 
with dj, a] by 

t fl ~ ^ (AS 

From a vacuum vector |0j) annihilated by a; one can construct a Fock space generated by 
the orthonormal set 



a i 

Ui ) = "Ti^rrl *)' Ui = '"' ' k ~ 1 > ( 5 ) 
where [n] 9 stands for the g-deformation of integer numbers 

and the factorial is defined by [0] 9 ! = 1, [n] q \ = — l] q ■ ■ ■ [l] q . In this representation 

one can readily express the action of the basic operators 

Ni\ni) = riilrii), 



a l\ n i) = \J\ n i + l ]q\ n i + 1), ( 7 ) 

a*K) = y^NgK - 1)- 

Being q = e 1 ^ a rational phase one has 

H, = s -^§ (8) 

so that [k] q = [0] q = and the Fock space ~ C fc is finite dimensional, with = 
^ajj = 0. Due to the symmetry [n] q — [k — n] q , one has [n] g ![A; — 1 — n] q \ = [k — l] q \ for 
any n = 0, • • • , — 1. The relations ([?]) are conveniently seen as finite dimensional matrix 
analogues of the usual harmonic oscillator (Heisenberg-Weyl) algebra, with g-deformed 
integers [n] q instead of integers n. From the deformed algebra (j2]), notice that one recovers 
the usual commutation relations 

[Ni, at] = -Oi, [N t ,al] = al (9) 

reflecting that Ni is indeed the number operator for the i th mode. It is simple to recover 
g-commutation relations from the matrix representation, by noting that [n + l] q = 
q[n] q + q~ n . 
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From the above facts, the fc-nilpotent g-bosons describe a simple system exhibiting ex- 
clusion statistics. They may be seen as hard core particles, with k regulating the core 
hardness: they combine properties of bosons, but carry in their very formulation a max- 
imum occupation constraint. It is apparent that in the limit k — > oo (i.e. q — > 1) one 
recovers standard bosons. However, the case k = 2 does not describe fermions (see (F2])). 

Regarding the commutation properties for different modes, we follow the criteria that 
g-bosonic operators corresponding to different degrees of freedom commute^—, 

[a h aj] = [a{, a]) = [a h a]] = for i ^ j. (10) 

This bosonic behaviour is set even for k = 2, another departure from fermions in our 
treatment. The Fock space of the system is thus simply "H = ®iHi- In comparison 
with standard bosons, recovered in the limit q — > 1, one must stress that a unitary 
transformation U (m) of g-bosons does not render g-bosonic modes^. 



III. COHERENT STATES - SINGLE PARTICLE CASE 

We discuss in this section one single mode i. Coherent states in Hi may be defined as 
eigenvectors of af, however, one readily notes that the only eigenvector of in the finite 
dimensional space Hi is the vacuum, as it happens for fermionic operators. The way 
out is to enlarge the Hilbert space by allowing linear combinations with coefficients that 
go beyond the complex numbers. One is lead to introduce fc-nilpotent para-Grassmann 
numbers, in the same way that Grassmann numbers are needed to deal with fermionic 
coherent states^. Consider an indeterminate 9i and a formal vector 

fc-i 

10;) = ^^k), (11) 

where a Ui are complex coefficients (we introduce the notation "| )"- to distinguish this 
expression from proper vectors in Hi), and evaluate 

k-2 

Oi\0i) = Oi an+iO?y/[ni + l} q \ni). (12) 

7li=0 

The conditions for having an eigenvector are 

a nt+1 = a n J J [rii + \} q (13) 
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and 

ef = o, (14) 

imposing the /c-nilpotency condition. One then gets 

k-1 , s 

i^) = E( a o/>/w)«r i in*> ( 15 ) 

satisfying 

a i \e l ) = e i \e l ). (16) 

Next, one introduces formal dual vectors by conjugation: consider 9{ an indeterminate 
conjugate to 0j and dual vectors (nj| in %* to define 

fe-l / X 

{0i\ = ^(m\[^/y/w)^ (17) 

so that 

(^|a| = [9 t \9 t . (18) 
The action of (0j| on |0j) gives a polynomial in 0j, 0, 

5«i /TO 

which is not a real number and cannot be normalized; we adopt the convention a — 1- 
We remark that, up to this stage, the construction includes standard Grassmann numbers 
for k — 2. 



Commutation relations and conjugation 

Before prescribing an iterated integration rule over 0j and 0j (see below), one needs a 
commutation relation 9-fii = a9i9i to be able to re-order general monomials. One usual 
criteria is to ask for nilpotency of linear combinations x% = P@i + oQ% with complex 
coefficients p, <r-, leading as the simplest solution for a a primitive complex root of one of 
order k, for instance a = e l ^~ = q 2 , 

0.0. = e^OiOt. (20) 

This is the usual choice in the literatures^^. Notice that only for k = 2 one finds that Xii 
Xi satisfy the same commutation relations as 0j , 0j (anti-commuting fermionic case, with 
real a = —1). For k > 3 Xi-, Xi wm n °t have the same commutation properties (1201) that 0j, 
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9i have. Then, in contrast with Grassmann numbers, a linear change of para-Grassmann 
generators with complex coefficients cannot preserve both k-nilpotency and commutation 
rules (l20]y ^. Moreover, the relation f l20|) has a serious drawback for k > 3: it does not 
support the usual conjugation of products (with the property (99') = 9' 9)&&. This has 
led some authors to avoid the use of conjugation^ or to adopt a non-standard conjugation 
rule for products^. 

We find no case in enforcing nilpotency under complex linear transformations while com- 
mutation rules of the resulting combinations are drastically different from that of the 
original ones. In this work we consider instead the relation 

9 l 9 l = a9A (21) 

with a G l£, with the advantage of supporting standard conjugation.— 

We remark again that with the usual commutation relations (}2T)j) one cannot make a 
linear transformation of para-Grassmann generators into para-Grassmann generators, in 
the sense of preserving nilpotency and commutation relations. In the multi-particle case, 
this excludes the use of Fourier transformations or any linear change of basis and makes 
it impossible to map interacting modes into decoupled ones, even for quadratic Hamilto- 
nians. Our choice (121 j) is neither better or worse in this sense, but allows for a consistent 
definition of conjugation and para-Grassmann symbolic calculus. Moreover, it leads to 
noticeable simplifications in the applications. 

Algebraic structure 

The construction discussed above contains the following algebraic structures: 

First an algebra C[0j, 6i\f (6^, 9f, 9 i 9 i —a9 i 9 i ) (quotient of the non-commutative free algebra 
of polynomials in 9i, 9i with the two-sided ideal generated by 9f, 9\ and 9fii—a9i9i ), which 
is a vector space of dimension k 2 over the field C with a closed product of vectors. This 
will be called^ the complex para-Grassmann algebra PGk, a , characterized by a nilpotency 
order k > 2 and a real commutation coefficient a (PG2-1 is the standard Grassmann 
algebra). Notice that it is also a ring, with the sum and product of polynomials. 

Second, the free module of the orthonormal set in %i over the ring PGk^- This is 

an extension of the Hilbert space Hi, the linear span of a basis over coefficients (numbers) 
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in PGk, a that are more general than complex numbers, and will be called 

= \ \ v ) = X 7"*! 71 *) such that G PG kta \ . (22) 

We will not distinguish left or right multiplication of PGk tU} numbers with vectors, so ZQ 
is technically a bimodule. 

As it is usual in the fermionic case, one can use functional analysis language calling 8i, 
§i para-Grassmann variables and writing the PGk, a algebra elements as functions of such 
variables 

/mo =£;/«»' w', (23) 

n,n' 

where the conditions Of = Q\ = Q$i — a0i8i = ensure that an arbitrary function is 
represented by k 2 complex coefficients (the form in (|23j) may be called an expansion 
of / in the anti-Wick ordered basis of PGk, a )- These functions are called holomorphic 
(antiholomorphic) when only powers of 9{ (9i) are present. Conjugation in PGk, a is then 
written as 

r(Mi) = £iw^ (24) 

n,n' 

where f nm stands for complex conjugation. The *-algebra property (fg)* = g*f* is 
fulfilled*. 

A sesquilinear form is naturally defined on /Q by extension of the inner product in 

given \v) = J2t~=o v m{ d h 6~i)\ n i) and \v) = Ej~i Vm(0i, we define 

fe-i 

{v\v) = Y,^MA)^MA)- (25) 

rii=0 

This is not an inner product (positivity does not make any sense), but is useful to write 
the projections of elements in /Q onto the basis §5§. According to previous notation, 

(m\u) = Vm(0i,9~i), 
(v\rk) = (rk\vT = %tM)- (26) 

In particular, 

{mm = (l/y/w) e?. (27) 
Integration 

Following Berezin's seminal work on Grassmann integration, one defines a linear form 
on PGk, a with integral-like properties. This, and a whole proposal for para-Grassmann 
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integral and differential calculus, has been done before using the commutation relation 
(r20j) 3 ' -i^. We will not pursue here such a complete program, that would presumably be 
simpler for PGk, a ~] we just quote the basic Berezin-like integration rules for anti-Wick 
ordered basis elements 

J e^' de t = e™N5 n ,^ (28) 

where M is a positive normalization constant, and we stress that 6i, 9i act as independent 
variables under integration. Then 

J d9i Bffi dh = Af 2 6 n ,k-iS n >,k-i (29) 

can be seen as a double iterated integral. The order of the factors and differentials must 
be cast as it is in ( |29|) before using the recipe. In what follows we set for convenience^ 

A/" = yj[k-l] q l (30) 



Completeness 



The completeness of the coherent states construction is expressed as a "resolution of 
unity" in One has to make sense of 

J d6i 10^(0*, 00(0*1 d6i — I\ Hi -> H h (31) 

where /ii(0j,0O may be seen as a measure weight for the integral (because of ordering 
issues it is important to set a position for the weight factor; we write it for convenience 
in the middle). 

Given any two vectors \u) = Y^=o u n.i\ n i) ^ \ v ) = Yln~=o v ^\ n i) in we then ask ^(0,, 0~O 
to fulfill 

(v\ f d9 i \9 i )^ i (9 u 9 i )(9 i \d9 i \u) = (v\u) (32) 



which amounts to 

QT ^ 



Mi -7=^(^, 9 i )^=d9 i = 8 ni< . (33) 

Writing jU»(0i, #0 in the general anti-Wick form (|23p . the expression under integration is 
anti-Wick ordered and may be solved with the rules f j29|) . The weight factor must contain 
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only terms with equal powers of 0, and 9i so as to produce non-vanishing results only for 
rii = n[. One gets 

fe— 1 

^(MO = £ ( 34 ) 

P =o Viv- 
as the unique (anti-Wick ordered) kernel making sense of (j3Tj) . 

Equation (13T]) singles out the auxiliary role of para-Grassman numbers and module /Q 
vectors in our construction: when computing matrix elements in Hi, module vectors 
are projected onto "H, and the expression leads to compute a complex valued integral 
on PGk, a - Namely, para-Grassmann numbers are "integrated out" to recover results in 
the g-boson Fock space, as it happens with standard Grassmann numbers in fermionic 
theories. 



Anti-normal order prescription 

Once we set 9i) as a measure weight, the use of the identity resolution (THTj) may lead 
us to integrals where the factors of 9i, 9% are not anti-Wick ordered. Following^ we define 
a linear anti- normal order prescription : :, moving in each term under : : all 9i factors 
to the left and all 8i factors to the right, without using commutation rules. 

This prescription is useful in several situations. First, the weight factor in ( 134)) can be 
written as 

^i(Mi) =: ef 4 : (35) 
where the g-deformed exponential is defined, as it is usual in g-deformed algebras, by 

< = E^'- < 36 > 

Second, one can define non-ambiguous Toeplitz operators from a PGfc iCr -valued symbol: 
given a function (j)(9i, 9j) one considers homomorphisms in Hi of the form 

T mA) = JdOi-. |0OMMiMMi)(0il : Mi ■ Hi Hi. (37) 

In this notation, these are generically called anti- Wick or contravariant operators^ with 
symbol 4>(9{, 9i) and are intimately related to Toeplitz operators^. For bosonic coherent 
states | z) associated to Lie algebras^ such operators are "diagonal" in the coherent states 
basis^; once applied to Fock space states, and projected onto the coherent state basis, they 
properly become Toeplitz operators mapping holomorphic square integrable functions on 
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a Kahler manifold onto themselves, through a Bargmann projection. They implement the 
Berezin-Toepliz (or coherent states) quantization of the classical function (p(z,z). 

The same structure may be realized here; indeed, ordering ambiguities and the conflict 
between para-Grassmann conjugation and commutation relations in (1201) prevent a con- 
sistent Berezin-Toepliz quantization of a para-Grassmann algebra. The ordering problem 
was solved in^ by introducing an anti-Wick ordering prescription, but still with commuta- 
tion relations as in (120]) requiring a non *-algebra conjugation. More recently^, the consid- 
eration of /c-nilpotent para-Grassmann algebras with independent (real) a-commutation 
relations as in ( l2Tj) allowed to construct a well defined reproducing kernel (expressing the 
Bargmann projection and therefore the resolution of the identity) and Toeplitz operators. 
It is noticeable that the g-boson operators and a] can be written as the Berezin- Toeplitz 
quantization of the simple symbols 9{ and 9i, respectively^. These recent papers are the 
basis for our present approach. 

A sesquilinear form is naturally defined in the para-Grassmann algebra^ as 



with f(9i,9i), g(9i,9i) in PGk, a - in contrast to (12 5p . this definition does provide an inner 
product in PGk, a - 

We remark that under any order prescription, the commutation rules for para-Grassmann 
variables play no further role; the /c-nilpotency, defining [n] q , and the measure weight 
/j,(9i,9i), setting orthonormality, are the key ingredients of the present construction. In 
particular, one can manipulate the use of (1ST]) in a clean way under the anti-normal order 
prescription. In Section |V] we develop simple trace formulae for operators acting on Hi. 



For a system with m independent degrees of freedom, the coherent states are the direct 
product of single mode coherent states. The key point in this Section is that, in our 
scheme, handling multi-particle coherent states presents no further complications. 

Formally, we first introduce m complex para-Grassmann variables 9i, ■ ■ ■ ,9 m . As the 
different mode operators commute (see (II Op ) and we do not require nilpotency of linear 
combinations of para-Grassmann variables, for i ^ j we set 




d9 i : r(9 i ,9 i ) f i(9 i ,9 i )g(9 i ,9 i ) : d9 t 



(38) 



IV. COHERENT STATES - MULTI-PARTICLE STATES 



9i9j — 




9{9j 



9j9i. 



(39) 
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We then define the m-mode para-Grassmann algebra P™ Q = 
C[0i,0i,--- ,9 m ,9 m ]/{9f,9f,9i9i - u9i9i,9i9j - 9j9i) as the quotient of the free al- 
gebra of polynomials in m complex indeterminates with the ideal expressing /c-nilpotency 
and all commutation relations, and the direct product of modules K, m — K,\ ® • • • ® K. m . 
Denoting 9 = {9i, • ■ ■ , 9 m }, we define coherent states as 

\9) = \9 1 )®---®\9 m ). (40) 



The elements of PG™ a are functions of several para-Grassman variables, that have unique 
coefficients when written in anti-normal order 

f(0,§)= /whC-^-? 1 , (41) 

{n}{n>} 

where {n} = {rii, ■ ■ ■ ,n m } is summed over each rii — 0, ■ ■ • , k — 1. Strictly speaking, 
the order is already set when each 9i is on the left of the corresponding 9f, we annotate 
a complete order of para-Grassmann variables, with decreasing indices, for convenience 
in solving integrals. For several variables we use for the anti-normal order prescription 
the same notation : : as before, moving in each term under : : all 9i factors to the left 
and all 9i factors to the right, without using commutation rules, and ordering commuting 
variables in decreasing index order just for convenience. 

Integration is defined iteratively. For the function in fj4T]) . the integral reads 



f d9 f(8, 9)d9= hn}{n'} fd9 1 --- d9 rn 9 n ™ ■ ■ ■ 9^ ! 9^ • • ■ 9>fd9 x ■ ■ ■ d9 m = 

{n}{n>} 

hn}{n>} J d9 l 9 n ^--- J d9 m 9 n ™ 1 9<d9 m --- J 9fd9 x (42) 



W{n'} 



providing a non vanishing result, from ( 12"5|) . only from the term with n\ = n' x 



n m — n 'm — k 1. 



Resolution of the identity in Tt — Hi <S> ■ ■ ■ <8> Hm is readily written as 

J d9 \6)n(8, 9){9\ d9 = I : H -> U, (43) 

with measure weight 

H(9, 9) =: e e q 1§1 e 9 ™ Sm :=: e 9 / 1 ■ ■ ■ e e f m : (44) 

We remark that the multi-particle version is a simple generalization of the one particle re- 
sults. This is due to the normal order prescription, partially taken frorn^. Indeed, we have 
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changed the g-commutation rules proposed by these authors for different para-Grassmann 
variables, which requires an extra order prescription and generates a conjugation problem. 

It is simple to generalize the results in& to multi-particle states. A sesquilinear form in 
PG% a is defined by 

(f(9,9),g(9,9)) = J ' dB : r(9,9)^(9,9)g(9,9) : d9. (45) 

An anti-Wick operator 

T m S) = fdd : \9M9,9)<j>(9,9)(9\ : dO (46) 

can be projected onto coherent states defining a Toeplitz operator-; in this way, creation 
and annihilation operators act on holomorphic functions^ by 



a-i = Tl 



1 ■ 



4 = T S , (47) 

Our objective here is to apply the above formalism in the construction of coherent states 
trace formulae. 



V. TRACE FORMULAE 

Given an operator A : Ti — > H, its trace can be written as an integral over coherent states 
in much the standard way. In the one particle case, say for a mode i, 

k-l 



";=0 



Tr(A) = ^2(rn\A\ni) 

ni=0 
k-l r 

= d9 i (n i \A\9 i )ii i {9 i ,9 i ){9 i \n i ) d9, 

fc-i 

d9i : Hi{9iA)Y^{8i\ni) H A%) : d9 { 

71=0 

= J d9, : lHiPiA){Pi\A\di) : d9 it (48) 

where we used the identity (13T1) and a reordering of factors under the anti-normal order 
prescription. 

In the multi-particle case we can operate the same way, using the identity (|43p and 
commutation relations (I39p . We start writing the trace in the canonical basis |{^}} = 
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\n x ) <g> • ■ • ® \n m ), 



Tr{A) = £<{n}|.A|{n}> 

{n} 



J d9 {{n}\A\9)fi(9,9)(9\{n}) d9 

{rii} 

dO : n(9,9)J2(9\{n})({n}\A\9) : S 
{„} 

= j d6 : n{0,6){9\A\0) : S. (49) 

As said before, no further complications arise in handling multiparticle states in terms of 
independent para-Grassmann variables. 



VI. APPLICATIONS: THERMODYNAMICS IN SIMPLE EXAMPLES 

We are interested in computing the thermodynamical partition function for a system of 
nilpotent g-bosons with Hamiltonian H at temperature fc#T = 1/(3. We thus need to 
evaluate coherent state matrix elements (9\e~^ H \9), a task that provides closed results 
only for some simple Hamiltonians. 

A. One q-boson in a thermal bath 

Consider one g-boson oscillator (q = e™l k ) with Hamiltonian 

Hx = eN h (50) 

having spectrum e ni = nie, n\ = 0, • • • , k — 1 (notice that here n\ counts the excitations in 
the one particle spectrum, not particle number). We compute the canonical (one particle) 
partition function in a thermal bath, 

Z 1 (P) = Tr(e~^). (51) 

In order to compute (9i\e~ l3Hl \9±) it is convenient to expand the coherent states using 
([15]), getting 

nninni 

(0i|e-^|*i) = EtHT^" 1 - ( 52 ) 

ni=0 

The trace formula (J3HD is easily integrated using the rules (129"j) providing 

k-i _ kj3e 

Z 1 (/3)=^e-^ = T -^ r . (53) 

ni=0 
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Figure 1: (colour on line) Mean energy of a q-boson in a thermal bath, for different nilpotency 
orders k (in arbitrary units). Notice that at finite k the mean energy saturates at high enough 
temperatures, where each of the energy levels is found with equal probability. In contrast, for 
standard bosons (k oo) the mean energy grows linearly with temperature. 
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Figure 2: (colour on line) Specific heat for different nilpotency orders k (in arbitrary units). 
Notice that at finite k the specific heat has a maximum and decays to zero when the mean 
energy saturates. For standard bosons this does not occur (no saturation is possible). 

This is of course the trace result straightforwardly computed in the canonical basis ([5]). 
The corresponding mean energy reads 

dlogZt(P) ( 1 k 



EM 



df3 



3 /3e 



(54) 



and the specific heat 



c{p) = -p dEl 



dp 4^ t^sinh 2 {fie/2) sinh 2 (k/3e/2) 



(55) 



We show in Figures (pfl |2]) these functions for low values of k together with the limit case 
k — > oo, to make explicit that the fc-nilpotent behaviour interpolates between fermionic 
(k = 2) and bosonic (k — > oo) standard results. 

A similar analysis can be done for one g-boson oscillator with Hamiltonian 



H[ = ea\a\. 



(56) 
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which has spectrum e ni = [ni] g e, n\ = 0, • • • , k — 1. One obtains 

2 am am 

{6 x \e-^\e) = Y LJ-e-^, (57) 
so the partition function gives 

fc-i 

zm = J2 e~ Mm]q - ( 58 ) 



B. System of (/-bosons 

In setting a multi-particle system of nilpotent g-bosons one must take into account that 
a linear transformation (in particular the Fourier transformation) of g-boson annihilation 
or creation operators does not render modes with the same commutation relations. The 
same occurs with para-Grassmann variables in our approach (and any other in the litera- 
ture). We restrict to Hamiltonians in which different degrees of freedom do not interact. 
Speculatively, one can think of a system with a finite dimensional Hilbert space per degree 
of freedom, such as a spin S system in the presence of strong interactions, which after a 
suitable transformation leads to independent fc-nilpotent g-boson modes. 

Let us consider a system of m g-bosons aj, aj with Hamiltonian 

m 

H = J2 e i N r ( 59 ) 
i=i 

The grand partition function at finite temperature /c#T = 1//3 is given by 

Z{p) = Tr [e-m-^ ; (60) 
where N = Y^j=x ^Y? is the total number operator. One needs to compute 

m 

(9\e-^ H -^\9) = (6\Y[e-^-^\6) (61) 

3=1 

which simply factorizes to give 

m am am an m an m 

f = \ w [ni]! M 

The trace is computed according to 



/at, If* 

d9 :,,i0. 0)11(0, ( i' ^0,):<IO J] / <W , : MMi) E l^f ( : d0 3 

j=i j=i y nj L 

(63) 
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Figure 3: (colour on line) Average occupation number of levels e,- for a system of g-bosons at low 
temperature. For k = 2 the behaviour corresponds to fermions, while for k — > oo it corresponds 
to standard bosons. 

giving rise to 

™l_ e -fcj9(ej-M) 

^ = 11 ! _,-«.,-,) • < 64 > 
i=i 

The relevant quantity to compute here is the mean occupation number of levels ej, which 
reads 

n iW>ri = ePto-*)-l ~ eWi-ri-l - (65) 

For finite k this is well defined even for €j = \i (evitable singularity), while the limit 
k — > oo is finite only for 6j > \i (the correct behaviour for standard bosons). In Figure 
([3]), it can be seen that the mean occupation at ej = fi is (k — l)/2 (identical to fermions, 
for k = 2) and diverges for k — > oo (Bose-Einstein condensation). For k = 3 the present 
result for n(e) is markedly close to the distribution of Z 3 parafermions and of particles 
with g = 1/3 Haldane exclusion statistics illustrated 

We note again that the fc-nilpotent behaviour interpolates between fermionic and bosonic 
standard results. 



VII. CONCLUSIONS 



The construction of coherent states for g-commuting particles requires the introduction of 
para-Grassmann variables. In particular, when q is a complex rational primitive root of 
unity, q = e t7T ^ k , the required para-Grassmann variables are fc-nilpotent. Many attempts 
have been done towards a consistent formulation of nilpotent para-Grassmann calculus, 
leading to different difficulties in the multiparticle case. In consequence, no consensus has 
been reached yet in the proper characterization of such nilpotent variables. 
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We have traced back the source of difficulties, as well as figured forward the applicability, 
of para-Grassmann variables in systems of nilpotent g-bosons. In this work we present 
a construction, in line with recent proposals^, that incorporates on the one hand para- 
Grassmann commutation rules which are independent from the nilpotency order, and on 
the other hand a normal order prescription for the generalized Berezin integration. 

Our approach solves the conjugation problem for complex para-Grassmann variables and 
allows for a consistent symbolic para-Grassmann calculus. In particular it makes possible 
to handle in much the standard way a resolution of unity as a generalized Berezin integral 
of multi-particle coherent state projectors. This allows for simple trace formulae, which 
have been used here to study the thermodynamics of simple Hamiltonians; the distribu- 
tion of /c-nilpotent g-bosons in a multiparticle system turns out Fermi-like, with mean 
occupation per mode bounded by k. This exclusion statistics could find application in 
the study of the plethora of novel phases in strongly correlated systems, where different 
types of "novel" statistics have already shown up, mainly as a consequence of the strong 
interactions. Our formalism allows for a thermodynamical description, hence providing 
the tools to compare with experimental measurements to come. 
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